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In this paper, we study two-photon interference with the approach of photon quantum theory, 
with specific attention to the two- photon interference experiment carried out by Milena D'Angelo et 
al. (Phys.Rev.Lett87 : 013602, 2001). We find the theoretical result is accordance with experiment 
data. 

PACS numbers: 42.25.Dv, 42.82.Cr, 03.65.-w 
Keywords: Quantum theory; Two-photon interference 

1. Introduction 

Nonclassical interference is one of the most remarkable phenomena in quantum optics. In particular, 
it can be observed in experiments with spontaneous parametric down conversion (SPDC) [1], a nonlinear 
optical process in which high-energy pump photons are converted into pairs of low-energy photons (usually 
called signal and idler) inside a crystal with quadratic nonlinearity. It has been shown in many experiments 
that the quantum state of the signal- idler photon pair is entangled [2]. Many experiments have made 
use of SPDC to demonstrate fascinating topics in quantum optics, such as the test of Bells inequalities, 
quantum communication, quantum teleportation, etc. [3], and its possible applications include quantum 
communication, computation, and cryptography [4]. All these experiments basically belong to the same 
category: quantum interference. Two-photon interference is one of the pure quantum phenomena attributed 
to quantum correlations. In experiments on two-photon interference, each photon pair behaves like a quantum 
object called a "biphoton", whose effective energy (or frequency) is twice that of the original photons, and 
the interference fringe of the photon pair has half the period of a one-photon interference fringe. 

Two-photon interference is a powerful tool to study the fundamental problems of quantum theory. For 
example, the Einstein-Podolsky-Rosen problem [5] is believed to be resolvable by testing Bells inequality [6] 
and the Greenberger- Horne-Zeilinger theorem [7] in two-photon or multiphoton interference experiments. 
Two-photon interferometry also has broad applications in practical areas such as quantum cryptography 
[8], metrology [9], potentially in quantum computing [10], precision metrology, information processing and 
imaging Coincidence imaging, or ghost imaging [11, 12]. 

Recently it has been argued that classically correlated light might mimic some features of the entangled 
photon pairs in coincidence imaging setups. Notice that the possibility of simulating the two-photon imaging 
features of entangled states with classical sources was not ruled out by the authors of the original ghost imag- 
ing experiment [13]. Both the theoretical work of Abourraddy et.al. [14] and the experimental investigation 
of Bennink et.al. [15] stimulated a very interesting debate about the role of entanglement in two-photon 
coincidence imaging [16]. In this work, we study the two-photon interference with the approach of rela- 
tivistic quantum theory of photon. In the viewpoint of quantum theory, the light has the nature of wave, 
and it is described by wave function ip{f,t) for the photon of spin 1 . The absolute square \ip(f,t)\ 2 can be 
explained as the photon's probability density at the definite position. For light interference and diffraction, 
the interference and diffraction intensity / is directly proportional to \ip(r, t)\ 2 distributing on display screen, 
and the light wave functions can be divided into three areas. The first area is the incident area, where 
the photon wave function is a plane wave. The second area is the slit area, where the light wave function 
can be calculated by quantum wave equation of photon. The third area is the diffraction area, where the 
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FIG. 1: The single-slit geometry, a is the width and b is the length of the slit. 

light wave function can be calculated by the Kirchhoff 's law. For double-slit interference, we can obtain the 
total diffraction wave function by superposition the diffraction wave function of every slit. For two-photon 
double-slit interference, we calculate the total interference wave function ip s (f,t) = ciipi(r,t) + C2ip2{r,t) and 
ipi (f, t) — cz<p\ (r, t) + C402 (r, t) for the the signal and idler photon, and the detectors D\ and D 2 measure the 
interference intensities are directly proportional to \ip s (r,t)\ 2 and \ipi(f, t)\ 2 for the signal and idler photon, 
respectively. The intensity of coincidence measurement is directly proportional to \ip s (f,t) ■ ipi(r, t)\ 2 for 
two-photon double-slit interference. In the following, we shall calculate these wave functions, and compare 
the calculation result with the experiment. 



2. Quantum approach of photon single-slit diffraction 

In an infinite plane, we consider a single-slit, its width a and length b are shown in FIG. 1. The x axis 
is along the slit length and the axis is along the slit width a. In the following, we calculate the light wave 
function in the single-slit with relativistic wave equation. At time t, we suppose that the incident plane wave 
travels along the z axis. It is 



= E^ 



(pz-Et); 



= E* 



■ 0j -e * Et ej, 



(1) 



where ipoj — Aj ■ e« F , j — x,y, z and A is a constant vector. The time-dependent relativistic wave equation 
of light is [12] 



d - 

ih—ip(r, t) = chV x ip(r, t) + Vip(f, t), 



(2) 



where c is light velocity. From Eq. (2), we can find the light wave function ip(r,t) — > when V(r) — > oo. 
The potential energy of light in the single-slit is 



V(x,y,z) =0 0<a;<6,0<2/<a,0<z<c' 
= oo otherwise, 



(3) 



where d is the slit thickness. We can get the time-dependent relativistic wave equation in the slit (V(x, y, z) = 
0), it is 



d - 

ih—i>{r,t) = cfiV x ip(r,t), 
by derivation on Eq. (4) about the time t and multiplying ih~ both sides, we have 
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substituting Eq. (4) into (5), we have 

§p$(?,t) = -c 2 [V(V-^,t))-VV(r,t)], (6) 

where the formula V x V x B = V(V • B) - V 2 B. From Ref. [11], the photon wave function is ip(r,t) = 
L (E(r, t) + iacB(r, t)), we have 



V-^(f,i) = 0, (7) 

from Eq. (6) and (7), we have 

(|^- c 2 V 2 )V^) = 0. (8) 

The Eq. (8) is the same as the classical wave equation of light. Here, it is a quantum wave equation of 
light, since it is obtained from the relativistic wave equation (2), and it satisfied the new quantum boundary 
condition: when ip(f,t) — > 0, V(r) — > oo. It is different from the classic boundary condition. 

When the photon wave function -0(r, t) change with determinate frequency u>, the wave function of photon 
can be written as 

iP(r,t) = $(r)e- i " t , (9) 

substituting Eq. (9) into (8), we can get 

and the wave function satisfies boundary conditions 

$(0,y,z) = $(b,y,z) = 0, (11) 

ip(x, 0, z) = $(x, a, z) = 0. (12) 

The photon wave function vp(r) can be wrote 

ip{r) = ip x (r)e x + ip y (r)e y + tp z {r)e z 

= E V^e;, (13) 

j=x,y,z 

where j is x, y or z. Substituting Eq. (13) into (10), (11) and (12), we have the component equation 

^j{^,y,z) = i}j{b,y,z) = Q, (15) 



ipj(x, 0, z) — ipj(x, a, z) = 0. (16) 
The partial differential equation (14) can be solved by the method of separation of variable. By writing 

i> j (x,y,z)=X j (x)Y j (y)Z j (z). (17) 



From Eq. (14), (15), (16) and (17), we can get the general solution of Eq. (14) 
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since the wave functions are continuous at z = 0, we have 

ip (x,y,z;t) \ z=0 =^(x,y,z;t) \ z=0 , 

or, equivalently, 

ip 0j (x,y,z) \ z=0 = ipj(x,y,z) \ z=0 ■ (j =x,y,z) 
From Eq. (1), (18) and (20), we obtain the coefficient D mn j by fourier transform 



Dmnj = —r / A? sin— — sin d£dr) 
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= otherwise, (j = x, y, z) 

substituting Eq. (21) into (18), we have 
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wAx,y,z) = > -— — —sin sin 
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substituting Eq. (22) into (9) and (13), we can obtain the photon wave function in slit 

ip(x,y,z;t) = ^2 i>j( x ,y, z ,t)ej 
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We can consider the case of limit, i.e., the slit length b is infinity, and the Eq. (8) and (10) become 

d 2 - d 2 d 2 - 
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we can easily obtain the light wave function in the single-slit when b — > oo 
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3. The wave function of photon diffraction 

In the section 2, we have calculated the photon wave function in slit. In the following, we will calculate 
diffraction wave function, we can calculate the wave function in the diffraction area. From the slit wave 
function component ipj (f, t), we can calculate its diffraction wave function component $j (r, t) by Kirchhoff 's 
law. It can be calculated by the formula [17] 

^j(r,t) = ~ [ ^-lt ■ [vVj + (ik - h^r1>j]d8. (27) 
the total diffraction wave function is 

$(f,t) = Yl ^(r,t)ej, (28) 

j=x,y,z 

the diffraction area is shown in FIG. 2, where k = 2? is wave vector, s is the area of the single-slit, r the 
position of a point on the surface (z = c'), P is an arbitrary point in the diffraction area, and the n is a 
unit vector, which is normal to the surface of the single-slit. From FIG. 2, we have 

(29) 
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where k 2 = k'^. Substituting Eq. (22), (29) and (30) into (27), one can obtain 

p ikR r -> -v ~ . 
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Assume that the angle between k 2 and x axis (y axis) is ^ — a (| — /3), and a((3) is the angle between k 2 
and the surface of yz (xz), then we have 

k 2x = k sin a, fc 2a = £;sm/3, (32) 

it • &2 = kcos9, (33) 
where 6> is the angle between fc^ and z axis. Substituting Eq. (32) and (33) into (31) gives 
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Substituting Eq. (34) into (28), one can get 
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Eq. (35) is the total diffraction wave function in the diffraction area. From the wave function, we can obtain 
the diffraction intensity / on the display screen, we have 

J a \®(x,y,z;t)\ 2 . (36) 

4. Double-slit diffraction wave function of photon 

From Eq. (23), in the first slit, the photon wave function ipi(x, y, z; t) is 

v-^ ^ 164. (2n+l)TTX (2m+l)iry 
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FIG. 3: Double-slit geometry with a the single slit width, b the slit length and d the distance between the two slit. 

From FIG. 3, in the second slit, when we make the coordinate translations : 

x' = x 

y' = y-{a + d) 

z = z, (38) 
we can obtain the second slit photon wave function by the first slit light wave function. It is 

r . x , V- ^ lQA J {2n+l)irx . (2m + l)n[y - (a + d)} 
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With the Kirchhoff's law, similar as Eq. (35), we can get the light diffraction wave functions of the first 
and second slit, they are 



j—x.y.z m—0 n—0 



and 
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The total diffraction wave function for the double-slit is 

®(x,y,z;t) = c 1 S 1 (x,y,z;t) + c 2 S 2 (x,y,z;t), (42) 

where |ci| 2 + |c 2 | 2 = 1. From Eq. (42), we can obtain the counts C in the detectors D\ or D 2 is 

Ccx \$(x,y,z;t)\ 2 . (43) 

With the relativistic quantum theory of photon, we obtain the relation between diffraction intensity and the 
slit length, slit width, slit thickness, light wavelength and diffraction angle, which includes all interference 
and diffraction information. 
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In the experiment of two-photon interference [18], The 458 nm line of an Argon Ion laser is used to pump 
a 5mm BBO ((3 — BaB 2 0^) crystal, which is cut for degenerate collinear type-II phase matching to produce 
pairs of orthogonally polarized signal and idler photons. Each pair emerges from the crystal collinearly, with 
uj s « u>i s=s u> p /2, where uij (j — s,i,p) are the frequencies of the signal, idler and pump, respectively. The 
signal and idler are interfered by the same double-slit experiment device, and the interference-diffraction 
pattern of signal and idler photons are separated by the beam splitter PBS and are detected by the photon 
counting detectors D\ and D 2 , respectively. The output pulses of each detector are sent to a coincidence 
counting circuit for the signal-idler joint detection. The experiment setup is shown in FIG. 3 of Ref. [18]. 
Since the wavelength of signal and idler photons are equal, and the double-slit experiment device are same, 
their interference-diffraction wave function are same, i.e., Eq. (42). The counts of photon counting detectors 
D\ and D 2 are directly proportional to |$(x, y, z; t) | 2 , and the counts of photon coincidence counting detectors 
D is directly proportional to |$(x, y, z; t) ■ $(x, y, z; t)\ 2 . In the following, we shall compare the theory result 
with the experiment data. 

5. Numerical result 

The double-slit interference-diffraction experiment of two-photon had been reported by Milena D'Angelo in 
Ref. [18]. The experiment parameters are: the wavelength of signal and idler photons arc A s = Aj = 916nm, 
the width of each slit is a = 0.13mm, the distance between the two slits is d = OAmm. In theory calculation, 
we take the wavelength, the slit width and the distance between the two slits are same as experiment values. 
In calculation, the theory parameters are taken as: cl = 0.955, c2 = 0.298, A x = A y = A z = 0.896, the slit 
length b = 1.31 x 10~ 2 and the slit thickness c = 2.65 x 10 -5 . In solid curve is our theoretical calculation, and 
the dot curve is the experiment data [18]. From FIG. 4, we find that the theoretical result is in accordance 
with the experiment data, when the diffraction angle (3 is in the range of \ < 2(mrad). When the diffraction 
angle (3 is in the range of \(3\ > 2(mrad), the theoretical result has a small discrepancy with the experiment 
data. We think the experiment data should be measured accurately, and the theoretical calculation should 
be improved furtherly. 

6. Conclusion 

In conclusion, we study two-photon interference with the approach of photon quantum theory, and com- 
parison the theoretical result with the experimental data. We find that the calculation result is in accordance 
with the experiment data, when the diffraction angle j3 is in the range of |/3| < 2(mrad). When the diffraction 
angle (3 is in the range of \f3\ > 2(mrad), the theoretical result has a small discrepancy with the experiment 
data. We think the experiment data should be measured accurately, and the theoretical calculation should 
be improved furtherly. 
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FIG. 4: Comparison between theoretical prediction from (43) (solid line) and experimental data taken from [18] (circle 
point). 



